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Abstract 

f^ ■ Nonthermal radiation of a Kerr black hole is considered as tunneling of created particles 

CN . through an effective Dirac gap. In the leading seniiclassical approximation this approach 

bJO' is applicable to bosons as well. Our seniiclassical results for photons and gravitons do not 

^ ■ contradict those obtained previously. For neutrinos the result of our accurate quantum 

"^ , mechanical calculation is about two times larger than the previous one. 
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1 Introduction 

The amplification of an electromagnetic wave at tlie reflection from a rotating black hole, so 
called superradiation, was predicted by Zel'dovich [Q, and also by Misner [Q. Then it was 
studied in detail by Starobinsky and Churilov |^ for electromagnetic and gravitational waves 
(see also [§). It looks rather obvious that if the amplification of a wave at the reflection 
is possible, then its generation by a rotating black hole is possible as well. Indeed, direct 
calculation by Page [Q has demonstrated that the discussed, nonthermal radiation does exist, 
and not only for bosons, photons and gravitons, but for neutrinos as well. The last result looks 
rather mysterious since for fermions there is no superradiation. 

In the present work these processes are considered from another point of view: as tunneling 
of quanta being created through the Dirac gap. Certainly, this approach by itself can be valid for 
fermions only. It is clear however that in the leading semiclassical approximation the production 
of fermions and bosons is described by same, up to the statistical weight, relationships. 

Let us note that in the recent work by Calogeracos and Volovik [|| (we found out about 
it after the present paper had been written) an analogous mechanism was considered for the 
description of the friction experienced by a body rotating in superfiuid liquid at T = 0: the 
quantum tunneling of quasiparticles to the region where their energy in the rotating frame is 
negative. 

Our semiclassical results for photons and gravitons are in a reasonable qualitative agreement 
with the calculation [^. One cannot expect here a quantitative coincidence since for the most 
essential partial waves the action inside the barrier exceeds unity not so much if any. As to 
neutrinos, the imaginary part of the action for them, at the total angular momenta of real 
importance, is considerably smaller then unity. Therefore, for spin s = 1/2 we have performed 
a complete numerical quantum mechanical calculation of the effect. Here our result exceeds 
that presented in |Q by a factor close to two. Unfortunately, the lack of details in does not 
allow us to elucidate the origin of the discrepancy. 

2 Scalar field 

We will start with a problem of a methodological rather than direct physical interest, with the 
radiation of scalar massless particles by a rotating black hole. 

The semiclassical solution of the problem is started from the Hamilton- Jacobi equations for 
the motion of a massless particle in the Kerr field (see, for instance, |]^): 

dSr{r) y _ K^ [^r^ + a^)e-al,f 
Or J A ^ A2 ' ^ ^ 

dSem' , ( hV 

— = K. — \ ae smB . (2) 

Here Sr{r) and Se{6) are the radial and angular actions, respectively; 

A = r^ - r^r + a^; Vg = 2kM; 

a = J/M is the angular momentum of the black hole in the units of its mass M ; e is the particle 
energy; l^ is the projection of the particle angular momentum onto a (the velocity of light c is 
put to unity everywhere). 



As to the constant h? of the separation of variables, in the spherically symmetric limit 
a — i> oo it is equal to the particle angular momentum squared P, or to /(/ + 1) in quantum 
mechanics (the Planck constant h is also put to unity everywhere). The influence of the black 
hole rotation, i.e. of the finite a, upon k^ is taken into account by means of the perturbation 
theory apphed to equ. (||). The result is [Q] 
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Here and below we use the dimensionless variables lo = ekM, x = r/kM, 
us recall that in the semiclassical approximation the substitution 

/(/ + !)-.(/ + 1/2)2 

should be made. 

Besides, in the exact quantum mechanical problem, under the reduction of the radial wave 
equation to the canonical form 

R"+p'^{r)R = 0, 

in the expression for p'^{r) additional (as compared to the right-hand side of equ. (0)) non- 
classical terms arise, which should be, strictly speaking, included for / ^ 1. However, for the 
simplicity sake we neglect here and below these nonclassical corrections to p'^{r), which does 
not influence qualitatively the results obtained. 

The dependence of the classically inaccessible region, where the radial momentum squared 
p2 is negative, on the distance x is presented for different angular momenta in Figs. 1, 2. At the 
horizon the gap vanishes 0. For r ^ cxo the boundaries of the classically inaccessible region 
behave as ±(/ + l/2)/r. In other words, the centrifugal term for massless particles plays in a 
sense the role of the mass squared. Let us note that for / > 1 both branches of the equation 
p2(r) = fall down, but for I = 1 one branch near the horizon grows up, and the second one 
falls down. Thus, the radiation mechanism consists in tunneling, in going out of particles from 
the dashed region to the infinity. 

One should note the analogy between the emission of charged particles by a charged black 
hole and the effect discussed. In the first case the radiation is due to the Coulomb repulsion, and 
in the case considered here it is due to the repulsive interaction between the angular momenta 
of the particle and black hole . 

The action inside the barrier for the radial equ. (|I]) is: 
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the integral is taken between two turning points. For simplicity sake we confine for the time 
being to the case of an extremal black hole, a = kM. Let us note that due to a singular 
dependence of p on x, the action inside the barrier does not vanish at / > 1 even for the 
maximum energy a; = Izf^. So much the more, it stays finite at / = 1 (compare Figs. 1 and 2). 
The repulsive interaction is proportional to the projection Iz of the particle angular momen- 
tum and enters the tunneling probability through the exponent, but the barrier depends on 
the orbital angular momentum / itself. Therefore, clearly the main contribution to the effect 
will be given by particles with /^, close to /. The numerical calculation demonstrates that the 
contribution of the states with Iz 7^ / can be neglected at all. Besides, since the action inside 
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Figure 1: The energy gap for / = 1. 
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Figure 2: The energy gap for / = 2. 



the barrier decreases with the growing energy, the main contribution to the effect is given by 
the particles with energy close to maximum. 

Unfortunately, even for an extremal black hole an analytical calculation of the action inside 
the barrier cannot be carried out to the end. Therefore, to obtain a qualitative idea of the 
effect, we will use a simplified expression for k^: 
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(The results of a more accurate numerical calculation with expression (§) will be presented 
below.) In this approximation one can obtain a simple analytical formula for action inside the 
barrier for all angular momenta but / = 1. Let us assume that u = {1 — S) 1/2 with 6 <^1; just 
this range of energies contributes most of all into the radiation. Then, the turning points of 
interest to us, which are situated to the right of the horizon, are: 



3^1,2 = 1 + 



26 



2 ± Vl + 4// 



Now one finds easily 
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One can see from this equation that the term / (following P) in formula (|^) is quite essential even 
for large angular momenta: it generates the terms 4/1 in formulae (^ and (|^), thus enhancing 
15*1 for / ^ 1 by tt/a/S. Correspondingly, the transmission factor D = exp(— 215*1) gets about 
40 times smaller. Let us note that even the transition in k^ from /(/ + 1) to (/ + 1/2)^ makes 
the effect considerably smaller for / comparable to unity; but this suppression dies out for large 
angular momenta. 

It follows from formula (|^ that the action inside the barrier is large, it increases mono- 
tonically with the growth of /, starting with 15*1 = vr for 1 = 2. As to / = 1, one can see by 
comparing Fig. 1 with Fig. 2 that here the barrier is wider than for 1 = 2, and therefore the 
action should be larger. Indeed, the numerical calculation of the action inside the barrier IS"! 
with K^ given by formula (|^) confirms these estimates. Its results are presented in Table 1 
where for the comparison sake we present also the analytical estimates \San\ with formula (|^). 
By the way, this comparison demonstrates that the approximate analytical formula (H) works 



/ 


1 


2 


3 


s\ 


3.45 


3.15 


3.33 


t^an 




3.14 


3.34 



Table 1: Action inside the barrier for scalar particles. 



very well. The numbers presented in the table refer to an extremal black hole and maximum 
energy of emitted particles. It is clear however that the transition to nonextremal black holes. 



lower energies, and larger / will result in the growth of the action inside the barrier. Since here 

it proves to be always considerably larger than unity, the use of the semiclassical approximation 

inside the barrier is quite reasonable. 

Let us check now whether it applies to the left of the barrier. The corresponding condition 

is of the usual form: 

d 1 , . 

dxp{x) 

In other words, the minimum size of the initial wave packet to the left of the barrier should not 
exceed the distance from the horizon to the turning point. Near the horizon one can neglect in 
the expression for the momentum p{x) the term related to the centrifugal barrier, so that 

d 1 d {x-lf 
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One can see easily that for not so large /, which are of importance in our case, this expression is 
comparable to unity and condition (|^) does not hold. Nevertheless, despite of this circumstance 
and of the neglect of the nonclassical corrections to p^(r) mentioned above, the results of the 
semiclassical calculation, presented below, are correct qualitatively. 

Let us come back to the calculation of the radiation intensity. The radial current density of 
free particles in the energy interval de is at r -^ oo 

■ ( 7U - \^ ^^-P ^_\^ 2ndpr de 

Indeed, the initial summation here reduces to the integration over the azimuth angle of the 
vector I, which gives 27r, and to the account for the contributions of all projections Iz of the 
orbital angular momentum. As pointed out already, in our case it is sufficient to take into 
account only one of them, Iz = I- By means of the identity 

dpr = de, 



dpr 
we obtain in the result that the total flux of free particles at r — i> 00 is 
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One can easily see that in our problem the total flux of radiated particles differs from the 
last expression by the barrier penetration factor only. Thus, in our semiclassical approximation 
we obtain for the loss of mass by a black hole in the unit time the following expression: 
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Here the maximum energy of radiated quanta is 
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Th = km + y/k'^M'^ — a^ is the radius of the horizon of a Kerr black hole. The analogous 
expression for the loss of the angular momentum is 



^ = ^E /^exp(-2|5(e,/)|)&. 
^ 1=1 i 



(13) 



The results of the numerical calculation with formulae ( [TT| ) and (|IB|) of the loss by a black 
hole of its mass and angular momentum for different values of the rotation parameter a are 
presented in Table 2. We present here and below, for spinning particles, results of calculations 
only for sufficiently rapid rotation, a k, 1. The point is that with further decrease of a, not only 
the thermal radiation grows rapidly, but the effect discussed falls down even more rapidly. For 
smaller a this effect becomes much smaller than the thermal one, and thus its consideration 
there does not make much sense. 
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Table 2: Loss of mass (in the units of 10 ^ ttM'^) and angular momentum (in the units of 
10~^7rM), due to the radiation of scalar particles. 



As one can see from Table 2, the rate of loss of the angular momentum is higher, in the 
comparable units, than the rate of loss of the mass. In fact, it follows immediately from 
expression (|12|). Already from this expression one can see that even for the maximum possible 
energy the ratio of the corresponding numbers is 2:1. Real ratios are even larger. Hence an 
important conclusion follows: extremal black holes do not exist. Even if an extremal hole is 
formed somehow, in the process of radiation it looses the extremality immediately. 



3 Radiation of photons and gravitons 

The investigation of the radiation of real particles we start with the electromagnetic field. 
Photon has two modes of opposite parity: the so called electric mode, with / = j ± 1, and 
magnetic one, with / = j [|T0|. It follows from the duality invariance that the radiation intensities 
for these two modes are equal. Thus one can confine to the solution of the problem for the 
magnetic mode, and then just double the result. 

One can demonstrate that the situation with the gravitational waves is analogous. Again, 
there are two modes which, due to a special duality, contribute equally to the radiation, and 
for one of these modes I = j- 

Obviously, for a mode with / = j the radial equation in the semiclassical approximation is 
the same as for the scalar field, but with different value of k"^. It can be demonstrated also 
starting from the so called Teukolsky equation |[TT| (neglecting again nonclassical corrections 



to p^(r)). The corresponding eigenvalues of the angular equation for particles of spin s, found 
again in the perturbation theory, are [H]: 
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(14) 

We have included into this expression the term 1/4, necessary for correct semiclassical descrip- 
tion. Let us note that, as follows from the consideration of the helicity of a massless particle, 
the restriction j > s holds. Correspondingly, for a photon j > 1, for a graviton j > 2. As well 
as in the scalar case, the main contribution into the radiation is given by the states with the 
maximum projection of the angular momentum, jz = j. 

Let us discuss first of all whether the semiclassical approximation is applicable here. As 
to the situation to the left of the barrier, it does not differ qualitatively from the scalar case. 
The situation inside the barrier is different. As one can see from equation (|1^, the presence of 
spin makes k^ smaller, and correspondingly, makes smaller the centrifugal repulsion. In result, 
both the barrier and the action inside it decrease. This qualitative argument is confirmed by 
a numerical calculation of 151 for photons and gravitons with the maximum projection of the 
angular momentum j^ = j and maximum energy for the case of an extremal black hole (see 
Table 3). Therefore in the present case, one should expect that the accuracy of semiclassical 
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Table 3: Action inside the barrier for photons and gravitons. 



results is lower than in the scalar case. 

The semiclassical formulae for electromagnetic and gravitational radiation differ formally 
from the corresponding scalar ones (0) and (|13|) by extra factor 2 only, which reflects the 
existence of two modes. The results of this calculation are presented in Table 4. In it we 
indicate in brackets for comparison the results of the complete quantum mechanical calculation 
1^ which takes into account as well the thermal radiation. It is clear from Table 4 that even for 
a = 0.999, when the thermal radiation is negligibly small, our semiclassical calculation agrees 
with the complete one qualitatively only. It is only natural if one recalls that the semiclassical 
action in the present problem exceeds unity not so much, if any. This explanation is supported 
by the fact that for photon, where \S\ is considerably larger (see Table 3), the agreement of the 
semiclassical calculation with the complete one is considerably better. 
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Table 4: Loss of mass (in the units of 10 ^ ttM"^) and angular momentum (in the units of 
10~^ ttM), due to the radiation of photons and gravitons. 



4 Radiation of neutrino 

Let us consider at last the radiation by a rotating black hole of neutrinos, massless particles of 
spin 1/2. The wave function of a 2-component neutrino is written as (see, for instance P]): 

ip = exp{-iet + ij,4)) ( ^^ ^^ j . (15) 

It is essential that the wave equations for neutrino in the Kerr metrics allow for the separation 
of variables as well |Tl|. The radial equations in dimensionless variables are: 
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The angular equations are: 
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For K^ the same formula (|T^) takes place, but now of course with s = 1/2. As well as for 
bosons, it is sufficient practically to consider states with jz = j. 

It is essential that Ri corresponds at infinity, for x — *> oo, and at the horizon, ioi x ^ 1, 
to the wave running to the right, and R2 corresponds for x ^ 00 and for x — *> 1 to the wave 
running to the left. (For this classification, it is convenient to use the so called "tortoise" 
coordinate ^(x); for x — > 00 ^ :^ x — ^ +00, for x ^ 1 ^ ~ ln(x — 1) ^ —00.) It is quite 
natural that here the radial current density is 
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We are interested in the probability of penetrating the barrier for the state which is an 
outgoing wave at infinity. For a neutrino or antineutrino such a state has a fixed helicity. 



but has no definite parity. Meanwhile, the potential barrier depends in our problem, roughly 
speaking, on the orbital angular momentum, and therefore is much more transparent for the 
states of / = j — 1/2, than for the states of / = j ' + 1/2. (These states of given / have definite 
parity, and are superpositions of neutrino and antineutrino.) Moreover, at I = j — 1/2 for small 
j, which give the main contribution to the radiation, the action either has no imaginary part 
at all, or its imaginary part is small, so that our above approach is inapplicable. Therefore, we 
will solve numerically the exact problem of the neutrino radiation. 

Technically, it is convenient to find the reflection coefficient R in the problem of the neu- 
trino scattering off a black hole, and then use the obvious relationship for the transmission 
coefficient D: 

D = l-R. 

Here the expressions for the loss of mass and angular momentum by a black hole are: 
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The results, obtained by numerical solution of the system of radial equations ([TBI) , are presented 
in Table 5. In brackets we present the results of 0, which include the Hawking radiation 
contribution. For a black hole close to extremal one, at a = 0.99, where the thermal radiation 
is practically absent, our results are about twice as large as previous ones. 
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Table 5: Loss of mass (in the units of 10 ^ ttM^) and angular momentum (in the units of 

10"'^7rM), due to the radiation of neutrinos. 
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